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I. INTRODUCTION 

The linear-optical properties of metallic photonic crys- 
tal slabs (MPCSs) have recently attracted considerable 
attention because they can be viewed as a simple model 
system of two coupled oscillators: (i) a Lorentz oscillator 
electronic resonance couples to (ii) an electromagnetic 
resonance, (i) The electronic resonance comes about 
from charges which accumulate at the surface of the 
metal nanostructures when exposed to the electric field 
of the incident light. These charges induce a depolariza- 
tion field that can cither counteract or enhance the ex- 
ternal electric field, depending on the permittivity of the 
metal, hence depending on the frequency of light. The 
resulting resonance at the transition point is the well- 
known particle plasmon or Mie resonance^ (ii) The elec- 
tromagnetic resonance is the Bragg resonance of the pe- 
riodic arrangement with lattice constant a. Importantly, 
an appreciable coupling between these two oscillators re- 
quires an additional slab waveguide, e.g., underneath the 
metal nanoparticles. Therefore, the physics of metal- 
lic photonic crystal slabs is distinct from that of usual 
metallic gratings, which have been discussed extensively 
many years ago. 2 Tailoring the waveguide parameters al- 
lows one to control the coupling strength, since the cou- 
pling arises from the spatial overlap of the plasmon- and 
waveguide-mode fields. 

Two-dimensional MPCSs were first discussed in Ref.0 
employing gold nanoparticles on a dielectric waveguide. 



Later, 4 gold nanowires showed even more pronounced ef- 
fects. In the latter structures, the coupling of the inci- 
dent light to the particle plasmon resonance can conve- 
niently be switched on and off via the polarization. If 
the electric field vector is oriented perpendicular to the 
nanowires (TM polarization), a pronounced depolariza- 
tion field arises, giving rise to a strong optical resonance. 
In contrast, if the electric field vector is along the wire 
axis (TE polarization), the depolarization factor is zero 
and one rather gets a Drude-type response of the metal. 

More recently, nonlinear-optical experiments on 
MPCSs have been presented^ These were interpreted 
along the lines of similar experiments^ performed on 
metal nanoparticle ensembles on a substrate surface with- 
out a slab waveguide. Thus, our initial motivation was 
to continue with experiments along these lines and ob- 
tain additional information from nonlinear-optical exper- 
iments. 

This paper contains both theory and experiments and 
is organized as follows. In Sec.[n] we start by dis- 
cussing the nonlinear optics of metal nanoparticle en- 
sembles with inhomogeneously broadened plasmon reso- 
nances. We derive rigorous analytic results for excitation 
with S pulses and Lorentzian inhomogeneous broadening. 
Furthermore, we present numerical results for a Gaus- 
sian inhomogeneous broadening and finite-duration opti- 
cal pulses. For second- (SHG) and third-harmonic gener- 
ation (THG), we find that the interferometric nonlinear 
response exclusively depends on the total linear-optical 
linewidth of the particle plasmon, i.e., one cannot ob- 
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tain any information on the relative contributions of ho- 
mogeneous and inhomogeneous broadening, respectively. 
This finding is in striking disagreement with the claim 
of Ref. 6 that interferometric SHG together with linear- 
optical measurements can differentiate between homoge- 
neous and inhomogeneous broadening. Reference^ has 
been the basis of much if not most of the work that fol- 
lowed in this field AiSiSiiSiiiii^ The discrepancy is traced 
back to a simple technical mistake in that work^ where 
the authors have not properly differentiated between the 
contributions of second-harmonic generation on the one 
hand and optical rectification on the other hand. In- 
deed, optical rectification (OR), self-phase modulation 
(SPM), or four- wave mixing (FWM) would allow for such 
a differentiation. Next, in Sec. lIIII we discuss the linear- 
optical properties of two coupled oscillators (represent- 
ing the particle plasmon and waveguide resonance). In 
contrast to frequent belief, the optical response of two 
coupled classical damped Lorentzian oscillators does not 
correspond to that of two new effective Lorentzian os- 
cillators. Generally, one rather gets Fano-like lineshapes 
in the linear-optical spectra. In Sec. II VI we discuss the 
nonlinear-optical signals from two coupled oscillators. 
We show that signatures of interferometric THG depend 
on the source of nonlinearity. Our theoretical analysis 
thus yields additional insights compared to the discussion 
in Ref. 5. The parameters of the presented numerical cal- 
culations are chosen to allow for a direct comparison with 
our experimental results, which are presented in Sec.lVl 
Compared with previous work, our experiments are dis- 
tinct in two aspects, (a) and (b). (a) First, we use 5 fs op- 
tical pulses, which are within the range of the anticipated 
particle plasmon decay times of 0.7—9 fs»i£ Previous work 
used pulses of 13 fs duration and longeri^ (b) As usual 
in "time-resolved spectroscopy," indirect information on 
the temporal behavior is obtained by exciting the sam- 
ple with a pair of time-delayed pulses, e.g., in pump- 
probe or transient four- wave mixing experiments. It is 
known that additional information can often be obtained 
by spectrally resolving the probe beam or the diffracted 
beam. In analogy, one anticipates that spectrally resolv- 
ing the third-harmonic signal, generated by the sample, 
versus the time delay of two exciting pulses, gives addi- 
tional insight. Indeed, our experiments reveal that dif- 
ferent spectral components of the third-harmonic signal 
can exhibit substantially different temporal dynamics - 
information that would obviously not be available from a 
spectrally integrated experiment. The comparison of our 
experimental data with theory allows us to determine the 
dominant source of the underlying optical nonlinearity. 
Finally, we conclude in Sec. I VII 



II. NONLINEAR OPTICS OF ENSEMBLES OF 
LORENTZIAN OSCILLATORS 

To probe metal nanoparticles with diameters in the 10— 
200 nm range by linear- or nonlinear-optical techniques, 



one often averages over several thousands of these par- 
ticles in order to obtain an acceptable signal strength. 
Depending on the fabrication method (e.g., lithographic 
patterning^ Volmer- Weber growth^) , a distribution in 
particle size and shape results, leading to a distribution of 
plasmon-resonance frequencies^ Ensemble linear-optical 
experiments alone cannot distinguish this inhomogeneous 
contribution from the homogeneous linewidth (resulting 
from an expected plasmon decay time r of a few fem- 
toseconds). Therefore, e.g., a combination of linear and 
nonlinear methods has to be used to extract both homo- 
geneous and inhomogeneous contributions. Nevertheless, 
not all nonlinear methods allow for this determination. 

We begin by discussing analytic results for the limit of 
8 pulses and Lorentzian inhomogeneous broadening, and 
continue with numerical simulations for finite Gaussian 
inhomogeneous broadening and pulses of finite duration. 

A. Analytic calculations 

Following along the lines of Ref.0, we start by describ- 
ing the particle plasmon by an oscillating particle with 
charge q, mass to, and displacement x(t), driven by an 
electric field E{t) via 

x + 2-yx + nix + (£x 2 + r/x 3 + ...) = —E(t) . (1) 

TO 

For the interferometric experiments to be described, E(t) 
corresponds to a pair of copropagating pulses with time 
delay T. In linear optics, i.e., for £ = r\ = ... = 0, this 
leads to a Lorentz os cillator r esonance at the damped 
eigenfrequency fi = yf2jj — 7 2 with a half width at half 
maximum (HWHM) j = 1/T 2 = l/(2r), the homoge- 
neous linewidth. T% is the dephasing time. To first order 
in the laser electric field, the polarization PW oc x^\t) 
is given by 

^'(t)ocfl- 1 f dt' e-^'-*') sm[fl(t-t')]E(t'). (2) 
J — 00 

Upon excitation with resonant pulses, the Fourier trans- 
form of PW(t) contains frequency components around 
±0. To second order in the laser electric field, 
— €( x (t)) 2 is the driving term for the second-order dis- 
placement x^(t). Provided that this driving term is off- 
resonant with respect to f2, we obtain the second-order 
polarization P^ 2 ) (t) oc x^ (t) with 

a?W(t)«[xW(t)] a . (3) 

The Fourier transform of this expression contains fre- 
quency components around ±20, i.e., optical second- 
harmonic generation, and components around zero fre- 
quency, i.e., optical rectification. Furthermore, Ref.@ ar- 
gued (see also formulae in Refs. EollT^l that the signal S 
measured by a slow detector is given by the integral of 
the nonlinear intensity over time, i.e., 

/oc 
dt[P^\t)f. (4) 
-OO 
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It is crucial to note that this expression comprises both 
SHG and OR. This, however, is in contrast to what is ac- 
tually measured in a second-order interferometric auto- 
correlation (IAC) setup, where one selectively detects the 
SHG by means of a photomultiplier tube behind optical 
filters, which suppress contributions other than SHG. 17 
For reasons of simplicity and to allow for analytic re- 
sults, we first discuss excitation with a pair of 5 pulses, 



E(t) = E [S(t) + S(t - T)] 



(5) 



It is clear from the symmetry that the nonlinear signals 
only depend on \T\. Thus, we only consider T > in 
what follows. For a single homogeneously broadened os- 
cillator we obtain 

P (1 \t) oc ^- 1 {6(i)e- 7 *sin(fit) 

+Q(t - T) e" 7 ^^ nin[n(t - T)]} . (6) 

This leads to the second-order polarization 

pW(t) 

oc fr 2 (9(<)e- 27 *[l -cos(2fii)] 

+0(< - T) e- 27 (^ T ){l - cos[20(i - T)]} 



-2©(i - T) 



3 -7(2t-T) 



{cos(ftT) - cos[0(2t - T)]}) . 

(7) 



Let us now consider an inhomogeneously broadened 
ensemble of oscillators with fixed damping 7 and a 
Lorentzian distribution of eigenfrequencies f2 with dis- 
tribution function 



L/tt 



(8) 



which is centered around frequency f2. The HWHM 

of this inhomogeneous distribution is T. To work out 

the convolution, we approximate the prefactor l/Q 2 in 
2 

Eq. iJJJ by 1/57 . This approximation is justified in the 
limit r <C f2, which is usually well satisfied for lithograph- 
ically fabricated particles. These two steps together lead 
to 



P 



(2) 



inhom 



(*) 



dfi P (n)P i2) (t) 



oc +e(i)e~ 27t 

+e(i-T)e^*- T ) 

+2S(t - T) e-r(*-T)-rr cos( n T ) 

-e(t) e - 2 ( 7+r ) t cos(2nt) 

-Q(t - T) e -2(7+r)(t-T) cos[2 n ( i _ T )} 
-2Q(t - T) c -(7+r)(2*-T) cos p(2t - T)] . 

(9) 

The first three lines correspond to OR, the last three 
lines to SHG. Note that the latter solely depends on the 



total width 7 + T. In linear optics, the width of the in- 
homogeneous ensemble results from the convolution of a 
Lorentzian with homogeneous width 7 with a Lorentzian 
of inhomogeneous width L. This leads to a total width 
of the resonance in linear optics of 7 + T. Thus, both the 
linear response and the correctly calculated SHG depend 
in the very same manner on the homogeneous and in- 
homogeneous linewidth, and a distinction is strictly not 
possible. 

In contrast, the contribution from OR does not simply 
depend on 7 + L, potentially allowing for a distinction 
between homogeneous and inhomogeneous lincwidths. 
By erroneously including OR in the calculated interfer- 
ometric "SHG signal," one can seemingly separate the 
homogeneous and inhomogeneous contributions to the 
linewidth. 

We have performed an analogous calculation for 
the third-order nonlinear-optical response. For third- 
harmonic generation and T > 0, we find that the en- 
semble THG polarization is 



P 



(3) (t) 

THG, inhom v / 

oc -0(i)e- 3(7+r)t cos(3m) 

-3( 7 +r)(t-T) 



-0(t - T) B -»W+nK*-*) C0S [3r!(t - T)] 
-30(i - T) e -( 7 + r )( 3 *" T ) cos[n(3i - T)] 
-30(* - T) e -(7+r)(3t-2T) cos [n(3i - 2T)] . (10) 

The THG again only depends on 7 + L, and no informa- 
tion on the homogeneous linewidth 7 can be obtained. 
However, self-phase modulation would provide such in- 
formation. In a non-copropagating geometry, the lat- 
ter would give rise to a diffracted four-wave-mixing sig- 
nal. Corresponding calculations have been presented in 
Ref.El 

Broadly speaking, nonlinear-optical signals of the type 
oj + oj (SHG) or uj + oj + oj (THG), etc. do not allow 
one to distinguish between homogeneous and inhomoge- 
neous contributions to the linewidth, whereas signals of 
the type u - u) (OR) orw + w- w (SPM or FWM), 
etc. do allow for such distinction. The "-" sign in OR, 
SPM, FWM, etc., effectively reverses the time axis in 
analogy to phase conjugation. For example in FWM, the 
"-" sign leads to the well-known photon-echo response. 19 
At this point, a decay of the ensemble polarization due 
to inhomogeneous broadening (just interference) can be 
reversed, whereas damping due to homogeneous broad- 
ening (a dissipative process) cannot be reversed. 



B. Numerical calculations 

The presented analytical calculations for S pulses are 
appropriate if the (complex) laser electric field spectrum 
exhibits negligible variation on the scale of the homoge- 
neous linewidth 7. For longer pulses, we have performed 
numerical simulations. As described above, the correct 
way to calculate the SHG contribution is to spectrally 
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FIG. 1: Recalculated data of Fig. 2 of Ref.|fl The dashed line 
describes the simulated FWHM of the second-order interfero- 
metric autocorrelation (IAC) from a single resonant Lorentz 
oscillator (the particle plasmon) versus plasmon decay time 
r = T2/2. The solid lines represent ensembles of oscillators 
with eigenfrequencies following a Gaussian distribution, for 
which the width is determined by fixing the total extinction 
linewidth AA (squares, AA=70nm; diamonds, AA=80nm; 
triangles, AA=125nm). The autocorrelation width is indi- 
cated by the full symbols for considering only the contribu- 
tion of SHG. Corresponding results for (erroneously) includ- 
ing SHG and OR are shown by the open symbols. 



filter the second-order response of the oscillator ensem- 
ble. To obtain the final IAC signal as a function of the 
time delay T, the square modulus of this filtered second- 
order polarization has to be integrated with respect to 
frequency. To allow for a direct comparison with the re- 
sults of Ref.0, we also use sech 2 -shaped 15 fs pulses with 
a center wavelength of 780 nm, resonantly exciting the 
ensemble. The latter has a Gaussian distribution of res- 
onances and is discretized in steps of 1 nm. 

Figure^ shows the resulting full width at half max- 
imum (FWHM) of the interferometric autocorrelation 
as a function of the plasmon decay time r. The full 
symbols correspond to the correct calculation, whereas 
the open symbols erroneously comprise the OR contri- 
bution and qualitatively reproduce the results of Ref.@ 
(see their Fig. 2). For each of the curves in our Fig.^ 
the total linewidth of the linear-optical spectrum is fixed. 
The squares, diamonds, and triangles correspond to a 
fixed extinction linewidth (FWHM) of AA = 70, 80, and 
125 nm, respectively (see parameters of Fig. 2 of Ref.0). 
The dashed curve corresponds to a single (homogeneously 
broadened) oscillator for reference. The correct results 
and those including the OR contribution differ strongly - 
as in our analytical calculations. In particular, the slopes 
of the correct curves in our Fig.^ are very nearly zero 
(within typical experimental error bars of 1 fs), while the 
incorrect simulations have a small positive slope. Thus, 



using the correct curves one cannot infer the plasmon 
decay time from measured interferometric autocorrela- 
tions, whereas the incorrect curves erroneously suggest 
this possibility. 6 We conclude that, under inhomogeneous 
conditions, the homogeneous linewidth cannot be deter- 
mined by analyzing linewidths from linear optics and 
SHG (or THG) measurements. Consequently, we will 
refrain from making any quantitative statements about 
plasmon decay times from now on. 

We note in passing that the IAC acquires artificial 
"wings" 6 if the OR contribution is erroneously included. 
Indeed, such wings are visible in Fig. la of Ref.|6|. They 
disappear in the correct calculation (not shown). 

III. LINEAR OPTICS OF TWO COUPLED 
LORENTZIAN OSCILLATORS 

In this section, we discuss the linear-optical properties 
of two coupled Lorentz oscillators. As already mentioned 
in the introduction, this system can serve as a simple 
model for metallic photonic crystal slabs. The results of 
this Sec. can be compared with the linear-optical exper- 
iments (Sec. 01 and, furthermore, are the basis for our 
discussion of the nonlinear-optical properties in Sec. II VI 

Generalizing Eq. Q to two coupled, oscillating parti- 
cles of equal mass m leads to 

Xpi + 27pi±p 1 + f2p 1 :Epi + (NL)pi-f^2; W g = -^-E(t) , 

(11a) 

i ; wg + 27 wg ± wg + f^ g :E W g + (NL) wg -^iz;pi = -^-E(t) . 

m (Hb) 

Here, x p \(t) and £ wg (i) are the displacements represent- 
ing the plasmon and waveguide oscillations, respectively. 
The resonance frequencies, (homogeneous) half widths at 
half maximum, and oscillator strengths of the uncoupled 
system are denoted by Qy, jj, and qj (j = pi, wg), re- 
spectively. fi 2 represents the coupling strength between 
the oscillators. The nonlinear terms (denoted by NL) are 
discussed in Sec. lIVI and ignored here. 

In order to make the resulting formulas transparent, 
we immediately discuss a few parameters in terms of 
their experimentally relevant values. Since the uncoupled 
waveguide resonance is extremely sharpSS as compared 
to the plasmon width, we set the waveguide damping 
7wg = 0. In the following, we derive formulas for an ar- 
bitrary waveguide oscillator strength q wg ; however, most 
aspects can already be understood in the simpler case 
q W g = 0. For typical sample parameters, \q wg \ <C |<7pi|> 
i.e., the area under the extinction curve of the (uncou- 
pled) waveguide mode is much smaller than that of the 
plasmon. 

In the frequency domain, Eqs. I|ll|l can easily be solved 
analytically. For monochromatic excitation, i.e., for 
E(t) = £0 e _,wi + c.c, this leads to the first-order dis- 
placements Xj(t) = Xj(w)e~ wt + c.c. and the polar- 
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izations P^\uj) 



N qj x^\uj). N is the density of the 



oscillators. Note that in the case g wg = 0, only x& (u>) 



contributes to the polarization. The total linear polar- 
ization becomes 



N qU-u 2 + £l 2 wg ) + 2q pl q wg n 2 c + ql g (-u 2 - 2iw 7pl + fl 2 ) 



m (-uji -2iLu lpl + n 2 pl )(-Lj2 + ni g )-ni 

The linear susceptibility x^H^O = P^'fa) / (sqEo) and the absorption coefficient 



En 



a(w) = (^/c )Im[x (1) (^)] = 



CKpl 



[(« 



na)( w 2 



n 2 i 



04]2 + 47^(^2 _ 2 g) 2 



(12) 



(13) 



r 



immediately follow. Sq is the vacuum permittivity, cq the 
vacuum speed of light, and a p \ — Nq pl / (2meoCoj p i) the 
maximum absorption coefficient of the uncoupled plas- 
mon oscillation. 

Examples of absorption spectra are shown in Fig.^a) 
for q wg /q p i = and Fig.EJb) for q W g/q P i = +0.1. One 
obtains the anticipated anticrossing behavior. For g wg = 
0, absorption maxima appear at the spectral positions 

1/2 



We note that, e.g., for g wg — 0, the total absorption 
(|13[) can be rewritten as a sum of two "Lorentzians," 
but with strongly frequency- dependent dampings. In the 
time domain, these frequency-dependent dampings corre- 
spond to a non-Markovian (and non-exponential) decay. 
For f2 a < Q wg < Ob; one solution can be described by 



oscillator a with constant resonance frequency Q a 
frequency-dependent damping 



and 



n 2 u 

a.b 



n 2 ) 



12 ± (n 2 pl 



n 2 wg y/A + ni 



7a M 



(14) 

These positions coincide with the normal mode frequen- 
cies of the coupled, but undamped system. For small 
O c and for O p i = f2 wg , the corresponding Rabi splitting 



7pi" P i 
2a(u) 



i _2 '3 — 3 



"pi-Tpi" 2 



l/2\ 



LU < Q 



> o v 



is given by H, 2 /fl p i. Hence, 
considered as "resonant" if 



the two oscillators can be 

fipi - o wg | < nl/h with 



(15) 



and an analogous expression for the oscillator b. 



IV. 



NONLINEAR OPTICS OF TWO COUPLED 
LORENTZIAN OSCILLATORS 



= (Opi + fl wg )/2, and as "nonresonant" otherwise. In 
contrast to frequent belief, the lineshapes in Fig.[21do not 
correspond to the sum of two effective Lorentz oscillators. 
One rather gets a highly asymmetric, Fano-like lineshape. 
Usually a Fano resonance results from the coherent in- 
teraction of a discrete quantum mechanical state with a 
continuum of statesiSi* 2 ^ In our purely classical model, a 
single sharp oscillator coherently interacts with a strongly 
broadened second oscillator. The latter replaces the con- 
tinuum. One result of the Fano-like interaction is that 
one obtains zero absorption between the two absorption 
maxima. The position of this zero appears at the root 
of the numerator of (|13fl . i.e., at or near the spectral 

position of the (uncoupled) waveguide mode O wg . Intu- S eneral form is S iven b y the terms 



In this section, we discuss the nonlinear-optical prop- 
erties of two coupled Lorentz oscillators in terms of 
third-harmonic generation. We consider an inversion- 
symmetric medium, hence all second-order nonlinear 
terms in Eqs. (|llf> are zero. At first sight, one might only 
expect third-order nonlinear terms like (NL) p i cx x pl or 
(NL) wg oc x„ g in Eas. (|llfl . Mathematically, the most 



wg- 

itively, this minimum is a result of destructive interfer- 
ence, which effectively suppresses the response of the two 
absorption "channels," of which the polarizations have a 
phase difference near it. This phase difference will also 
be important in nonlinear optics (see Sec. II VI) . When q wg 
is changed from zero to a nonzero value, the positions of 
the absorption extrema shift slightly, and the two peaks 
exhibit different heights as an additional characteristic. 
A reduced absorption of the more waveguide- like channel 



%fc[Zplft>] [Zw B (t)] 



(16) 



appearing in the differential equation for Xj(t), respec- 
tively (j = pi, wg; k = 0, 1, 2, 3). Here, we are only inter- 
ested in THG, which is off-resonant. In a perturbational 
approach the THG contributions to the third-order dis- 
placements are given by 



results, e.g., in the case q v 
top curves in Fig.^b)]. 



r/q p \ > and fipi < n wg [see 



(*)« 



fc=0 



- k \x {1) 

rwg 



(*)]' 



(17) 
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FIG. 2: Optical absorption spectra according to Eq. 113H 
(solid lines). One observes an anticrossing behavior when 
varying the waveguide resonance frequency fi wg with respect 
to the fixed plasmon resonance frequency S7 p i. Note the highly 
asymmetric, Fano-like lineshape of the peaks. All curves are 
displayed on the same scale. <?wg/g P i equals in (a) and 
+0.1 in (b). Common parameters are hQ p i = 1.6 eV, M7 wg 
as denoted for each curve, 7 P i = l/(2r), r = 2.5 fs, and 
fi,fi c = 0.5 eV. The dashed lines represent the pure plasmonic 
absorption in the absence of coupling, i.e., for hfl c — 0. The 
gray areas shown for fi wg = 1.6 eV depict the square modulus 
of the waveguide amplitude, |:Zwj| 2 , each exhibiting a single 
peak. In (a), the square modulus of the plasmon amplitude, 



is roughly proportional to the corresponding absorp- 



tion spectrum. The vertical line is a guide to the eye. 



The eight parameters r/j^t can be reduced to four, i.e., 
Vk = J^jljVj.k with k = 0,1,2,3, because the optical 
polarization is given by the weighted sum of the displace- 
ments. This immediately leads to the following general 
form for the THG contribution to the third-order polar- 



P^(t) ^^^(f)] 3 - 1 ^^)] 1 . (18) 



fc=0 

We note in passing that this form is generally differ- 
ent from the ansatz P^(t) oc [P^(t)] 3 (in analogy to 
Ref.H), which leads to i lk oc (I) q^ k qi g . 

For the numerical computation of THG spectra, we 
start off in the time domain. E(t) is chosen^ to resem- 
ble the 5fs laser pulses of the experiments (5fs Gaus- 
sian pulses deliver qualitatively similar results for all 
conditions discussed below). Furthermore, we fix O p i = 
1.67 eV, O wg = 1.56 eV, O c = 0.66 eV, r = 1.06 fs, and 
9wg/?pi = +0.085. These parameters correspond to sam- 
ple A in Sec.E] which can be considered as "resonant" 
according to the definition given in Scc. lIIII Integration of 
Eqs. (|llfl yields the first-order displacements Xj(t) and, 
with l|18|) . the third-order polarization. The square mod- 
ulus of its filtered Fourier transform delivers the THG 
intensity spectrum. Spectra are calculated as a function 
of the spectrometer photon energy and the time delay 
between the two excitation pulses, T. 

We first discuss the case r]k oc 8k,o (Sk,i is the Kro- 
necker symbol). The corresponding data set is shown in 
Fig.jSfa). A cut at T — (not shown) reveals four broad 
but clearly distinct spectral peaks in the THG spectrum. 
The appearance of four peaks can easily be understood in 
the frequency domain, since the third-order polarization 
for this case is proportional to the twofold convolution 
of the displacement (to) with itself, this displacement 
containing two peaks [see, e.g., Fig.|2Ja)]. The relative 
weights of the four peaks can be estimated employing the 
time domain. Assuming S pulses, g wg = 0, and neglecting 
damping, the two effective oscillators (see previous sec- 
tion) have comparable amplitude, and the general form 
of the THG polarization is proportional to 



[cos(O a t) + cos(Obi)] 3 oc 



+ cos(30 a i) 
+3cos[(20 a 
+3cos[(O a -+ 
+ cos(30bi) 



h O b )i] 
2fi b )t] 



(19) 



This contains terms at three times the normal mode fre- 
quencies O a and Ob as well as spectral mixing products. 
The relative amplitudes 1:3:3:1 of the frequency compo- 
nents 30 a , 20 a + Ob, O a + 20b, and 30b lead to the 
intensity ratios 1:9:9:1. This means that the two central 
frequency components are more prominent, in agreement 
with the numerical findings in Fig.Ufa). 

The pronounced dips between the four spectral peaks 
are closely related to the Fano-like lineshapes discussed 
in Sec. IIIII In linear optics, the phase relation between 
the two effective oscillators (absorption "channels" ) leads 
to destructive interference, and hence to zero absorption 
in the dip. The same destructive interference is also re- 
sponsible for the deep dips in the THG spectra. 
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THG Intensity (lin.) 




Time Delay (fs) 

FIG. 3: (a) Optical THG intensity derived from the coupled 
nonlinear oscillators. The THG intensity is shown on a satu- 
rated gray scale, versus spectrometer photon energy and time 
delay T between the two excitation pulses. At T — 0, the 
THG spectrum exhibits four distinct peaks (the high-energy 
peak is amplified by a factor of 10 for the sake of clarity). 
The four peaks exhibit different temporal behaviors. Corre- 
sponding cuts at the spectral peak positions indicated by the 
white arrows in (a) are shown in (b). For better compari- 
son, the curves are normalized to the same maximum and are 
vertically displaced. Obviously, the first and fourth curves 
both have a smoothly decaying (upper) envelope, while only 
the second and third curves show an envelope resulting from 
a beating. The nonlinearity parameters used are rjk oc 5k, o- 
The other parameters are quoted in the text. Compare with 
the corresponding experiment (Fig.[H] below). 



THG Intensity (lin.) 




Time Delay (fs) 

FIG. 4: As Fig.|SJa), but for different nonlinearity parameters. 
The shown THG contributions result from a nonzero value of 
(a) rji, (b) r/2, and (c) 773 while keeping the other nonlinearity 
parameters zero. At T — 0, the number of spectral peaks 
is three, two, and one, respectively, (a), (b), and (c) are 
displayed on individual gray scales. 



The behavior of the THG intensity as a function of 
time delay T differs among the four spectral peaks. Cor- 
responding cuts at the spectral peak positions indicated 
by the white arrows in Fig.[3^a) are shown in Fig.^b). 
The curves exhibit the usual oscillations with the respec- 
tive fundamental and harmonic frequencies, enclosed in 
the (upper) envelope of interest. The first and fourth 
curves clearly show a smoothly decaying envelope for in- 
creasing \T\. In contrast, the envelopes of the central 
two curves (which are associated with the spectral mix- 
ing products) reveal a beating. In spectrally integrated 
measurements^ this distinction is not possible. 

So far, we have only discussed the case rjk oc 5k,o- 
Next, we calculate corresponding THG spectra for dif- 
ferent nonlinearity parameters (Fig.EJ. In each part of 
this figure, all nonlinear parameters are zero except for 
a single one. The parts (a), (b), and (c) result from a 
nonzero value of rj%, 772 , and 773 , showing three peaks, two 
peaks, and one peak, respectively. In the frequency do- 
main this can again be understood by the corresponding 
convolutions. Remember that x^\lu) contains two peaks 

for the values chosen here, whereas iwg (u>) only contains 
one peak (refer to gray areas in Fig.|2J. 

In general, all parameters r]k can have nonzero values 
simultaneously. When adding up the nonlinear contribu- 
tions to the polarization, interference can result in a THG 
intensity with amplified or suppressed spectral peaks and 
dips, spectrally shifted peak positions, or even with new 
peaks or dips which are not present at all in Figs.^Ja) 
and 2] We will not go into a detailed analysis. We only 
note that for 770 3> 771 > 772 = 773 — 0, the tendency is 
to suppress the high-energy peaks compared to the case 
rjo ^ and 771 = 772 = 773 = 0. 

The key feature of the calculations presented so far is 
that the THG spectra depend on the underlying source 
of the nonlinearity, i.e., they depend on which of the co- 
efficients rjk is nonzero. In other words, observing four, 
three, two, or just one peak in experimental THG spec- 
tra allows one to learn something about the system by 
comparison with theory. This, however, is only possible 
for a certain regime of coupling between the two oscilla- 
tors, which we shall refer to as the regime of "moderate 
coupling." Obviously, for very small coupling strengths, 
i.e., for small values of fi c , the four spectral peaks in the 
THG spectrum of the case rjk — 5k,o (discussed above) 
merge into a single peak. In the other limit, i.e., for large 
values of fl c , also x^ g (uj) exhibits several peaks (unlike 
the gray areas in Fig.|2Jl, which can, for example, lead 
to several spectral peaks in the THG spectrum for the 
case 77fc = 6k,a as well. By numerical calculations for the 
"resonant" case (i.e., O p i = Q wg — f2), for q wg = 0, and 
assuming 5 pulses, we can specify the regime of "moder- 
ate coupling" by the condition 0.15 < ^l 2 c /{2-yVL) < 1.35. 
Thus, if one wants to learn something from the compari- 
son of experiment and theory, the coupling parameter of 
a sample has to be tailored correspondingly. 
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FIG. 5: (a) Scheme showing gold nanowires on top of a set of 
dielectric layers forming the slab waveguide. Hf02 is used as a 
high-index material, while MgF2 serves as a low-index spacer. 
The polarization of the normally incident white light or laser 
pulses is perpendicular to the wires (TM polarization) for the 
experiments shown in Figs.|SJ|HJ and[5] Samples with different 
lattice constant a, wire width w, and spacer thickness d are 
investigated, (b) Scanning electron micrograph of the gold 
nanowires (light gray) on top of the waveguide (dark gray). 



V. EXPERIMENTS 

The system of interest, a metallic photonic crystal 
slab, is schematically shown in Fig.[5fa). The coupling 
strength fl 2 between the particle plasmon resonance and 
the Bragg resonance (waveguide mode) can conveniently 
be tailored by the spacer thickness d. It is clear that an 
increasing spacer thickness leads to decreasing coupling. 
We experimentally find that when choosing d — 30 nm, 
the samples are within the regime of "moderate coupling" 
defined in the previous section, with a normalized cou- 
pling strength of = 0.43 for sample A and 0.28 
for sample B (see below) . These two selected samples are 
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presented in the following as examples for the "resonant" 
and "nonresonant" cases, respectively (see definition in 

sec.rrm. 



A. Sample fabrication and linear-optical 
experiments 

First, the dielectric layers shown in Fig.OJa) are de- 
posited in a high-vacuum chamber at pressures around 
10~ 6 mbar via electron-beam evaporation. We use 
hafnium dioxide (Hf02) as the high-index material (n — 
1.95) forming the core of the slab waveguide between the 
quartz substrate (n = 1.46) and the magnesium fluoride 
spacer layer (MgF2, n = 1.38). These dielectrics have 
been chosen for their transparency in the total spectral 
range of interest as well as for minimum THG genera- 
tion from the dielectric layers (as we investigated in in- 
dependent experiments). The 5-nm-thin indium tin oxide 
layer (ITO, n = 1.9) is necessary to avoid charging effects 
in the electron-beam writing process. Next, a photore- 
sist layer is spun onto the sample, exposed by means of 
electron-beam lithography, and developed. Finally, a 15- 
nm-thin gold film is evaporated and the metal on the re- 
maining photoresist areas is lifted-off . Each of the result- 
ing gold nanowire arrays covers a total area of (60/im) 2 . 
The electron micrograph in Fig.^fb) shows an enlarged 
view of a typical sample, revealing the high quality of the 
resulting structures. Typically, we fabricate entire sets of 
arrays on one glass substrate. In such a set, e.g., the lat- 
tice constant a is varied from 500 to 650 nm in steps of 
25 nm, and the nominal wire width from around 120 to 
around 220 nm in steps of 20 nm. In this fashion, we fab- 
ricate and investigate a total of 42 nanowire arrays on 
each substrate. 

To connect to theory, the measured extinction spec- 
tra (negative logarithm of the intensity transmittance, 
referenced to the substrate without gold structures) for 
TM polarization and for normal incidence are compared 
with Eq. l(T3|) derived in Sec. lIIII (see Fig.EJ). In these 
experiments we use a white-light source focused with a 
numerical aperture (NA) of about 0.025. Using a yet 
smaller NA tends to make the extinction dip even more 
pronounced. A careful discussion of this aspect can be 
found in Ref.|3 We find a good qualitative agreement 
of our simple theoretical model and the experiments. 
From a least-squares fit of the theory to the experiment 
(see Fig.EJ we obtain all relevant parameters, leaving 
only the nonlinear coefficients rjk as free parameters for 
the nonlinear-optical experiments to come. The exper- 
imental parameters of sample A (B) are: a — 550 nm, 
w = 185 ± 5 nm (a = 625 nm, w = 195 ± 5nm). The 
fit parameters of sample A (B) are: hfl p \ — 1.67 eV, 
M7 wg = 1.56 eV, hn c = 0.66 eV, r = 1.06 fs, and 
q^e/qpi = +0.085 (M7 p i = 1.65 eV, 7jQ wg = 1.39 eV, 
Ml c = 0.54 eV, r = 0.97 fs, and g wg /<? P i = +0.049), 
where r = 1/(27). The additional fit parameter a p i, 
together with the coefficients r/k, determines the abso- 



O 

o 

X 

LU 



1.0 



0.5 



c 1.0 
o 

'o 

is 05 



(a) Aa 




(b) 





1.0 1.2 1.4 1.6 1.8 

Photon Energy (eV) 



2.0 



2.2 



FIG. 6: Extinction of the selected samples (a) A and (b) B. 
The thick curves show the extinction measured with a white- 
light source, referenced to the substrate without gold struc- 
tures. The thin lines are absorption spectra derived from the 
model of coupled Lorentz oscillators, Eq. 1131 . obtained by a 
nonlinear least-squares fit to the corresponding experimental 
data. The gray area in (b) depicts the electric field spectrum 
(square root of the measured intensity spectrum) of the laser 
pulses used for the nonlinear-optical experiments. 



lute strength of the THG signals. Sample A is resonant, 
i.e., h\fl p i - n wg \ = 0.11 eV < 0.27 eV = hfl 2 jh, while 
sample B is nonresonant, i.e., h\fl p \ — fl wg \ = 0.26 eV > 
0.19 eV = M7 2 /f2 ( see discussion in previous section). 



B. Nonlinear-optical experiments 

In our THG experiments, we use 5fs laser pulses de- 
rived from a laser system closely similar to the one de- 
scribed in Ref.l25l (81 MHz repetition frequency). The 
pulses are sent into a Michelson interferometer, which 
is actively stabilized by means of the "Pancharatnam 
screw."— The linearly polarized pulses emerging from the 
interferometer are focused onto the samples (normal in- 
cidence and TM polarization) by a spherical mirror with 
a focal length of / = 100 mm. To estimate the inten- 
sities in the spot and to determine the effective NA, we 
have measured the spot size and the Rayleigh length by a 
knife-edge method in the horizontal (vertical) direction. 
The determined spot radius of 12.5 (12.3) /im is signifi- 
cantly smaller than the size of the nanowire arrays and 
leads to a pulse intensity around / = 3.9 x 10 10 W/cm 2 
and a laser fluence of 200 /iJ/cm 2 used in the experiments 
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FIG. 7: Characterization of the 5fs laser pulses, (a) Second- 
order interferometric autocorrelation (SHG-IAC) measured 
directly (thick line) and calculated (thin line) from the mea- 
sured pulse spectrum under the assumption of a flat spectral 
phase, (b) Third-order interferometric autocorrelation (THG- 
IAC) measured directly (thick line) and calculated (thin line) 
in analogy to (a). 



described below (for an average power of 80 mW in front 
of the sample). We will argue later that this fluence is 
still within the third-order perturbational limit. The cor- 
responding Rayleigh lengths of 490 (560) /im lead to an 
effective NA of 0.025 (0.022). 

In Fig.0 we depict the characterization of the laser 
pulses. The thick curve in (a) shows the usual second- 
order autocorrelation obtained from a very thin )3- 
barium-borate SHG crystal. The thin line is the auto- 
correlation as calculated from the measured laser spec- 
trum [see gray area in Fig.[|Jb)] under the assumption 
of a spectrally flat phase. The good agreement indicates 
that the residual chirp on the 5 fs pulses is of minor im- 
portance. The thick curve in Fig.[7Jb) shows the third- 
order autocorrelation function measured via THG from 
the surface of a thick sapphire plate. The thin curve in 
(b) is the corresponding calculated response under the 
same assumptions as in (a). Again, the agreement is 
very good. Notably, the envelope of the THG signal has 
decayed by a factor of 4 for time delays of just two cy- 
cles of light. These curves in (b) can be considered as 
the apparatus function and have to be compared with 
the measurements to be discussed in what follows. The 
achieved ratio of 32:1 between the THG signal at zero 
time delay and large time delays, respectively, indicates 
good alignment of the interferometer. 

In the THG experiments, the emission from the sam- 



ples in the forward direction is collected by another 
spherical mirror (focal length / = 100 mm), spec- 
trally prefiltered by means of four fused-silica Brewster- 
angle prisms to suppress the overwhelming fundamen- 
tal laser light and spectrally resolved using a 0.5-m- 
focal-length grating spectrometer (with a grating blazed 
at 250 nm wavelength) connected to a uv-sensitive, 
back-illuminated, liquid-nitrogen-cooled charge-coupled- 
device camera. 

Figure^ shows a typical data set of sample A, contain- 
ing 600 individual spectra obtained in a total of about 
8min acquisition time. Here, the THG signal is plot- 
ted on a linear gray scale as a function of spectrometer 
photon energy and time delay. The exact same represen- 
tation has already been employed in the theory Sec. (see 
Figs. 13 and 0J|. Indeed, the linear-optical parameters of 
Figs.|3 and 0] correspond to those of sample A [com- 
pare linear spectra in Fig.|SJa)]. Obviously, the mea- 
sured nonlinear-optical spectra are much closer to those 
in Fig.|3|than to any of Fig.^J In particular, four peaks 
occur in the spectra at zero time delay. Also, the depen- 
dencies of the different spectral cuts versus time delay in 
Fig-EIb) closely resemble those in Fig.|3fb). Again, the 
envelopes of the first and fourth cuts show hardly any 
beating, whereas the envelopes of the second and third 
cuts reveal a pronounced beating behavior [note that 
the very weak fourth peak in Fig^SJa) spectrally over- 
laps with the wing of the third peak, resulting in a small 
residual beating] . This comparison allows us to conclude 
that the nonlinear model where only 770 is nonzero is the 
appropriate one. This means that the nonlinearity pre- 
dominantly originates from the particle plasmon - which 
is not a priori clear. This finding is consistent with the 
experimental finding that the nonlinear signal decreases 
by a factor of 19 (and the multi-peak features disappear) 
when going from the TM polarization used so far to the 
TE polarization. It is also consistent with the fact that 
the THG signals drop by a factor of about 20 when go- 
ing from the gold nanowire arrays to areas of the glass 
substrate where only the dielectric layers are present. 

It is important for our interpretation that the exper- 
iments are performed in the third-order perturbation 
regime - which is also assumed in the theoretical analy- 
sis. Higher-order contributions would obviously modify 
the ratio of 32:1 between the THG signal at zero time de- 
lay and that at large time delays. In the experiments, the 
ratio of 32:1 is reached within experimental uncertainty: 
From analyzing the upper (lower) envelope of spectrally 
integrated data like those shown in Fig.[S] but for time 
delays up to ±60 fs, we derive a ratio of 24:1 (35:1). The 
actual ratio - which refers to a comparison between zero 
and infinite time delay - must lie between these two ra- 
tios. 

In Fig.|51 we show the data set for the "nonresonant" 
sample B. As for sample A, four spectral peaks are vis- 
ible in the THG spectra. In contrast, however, the 
peaks in Fig.|^a) have rather different spectral widths, 
as expected from the fact that the two effective extinc- 
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FIG. 8: (a) Measured THG intensity from sample A. The 
THG intensity is shown on a saturated gray scale, versus spec- 
trometer photon energy and interferometric time delay. Near 
time delay T — 0, the THG spectrum exhibits four peaks (the 
weak high-energy peak is amplified by a factor of 40 to become 
visible). The four peaks exhibit a different temporal behavior. 
Corresponding cuts at the spectral peak positions indicated 
by the white arrows in (a) are shown in (b). For better com- 
parison, the curves are normalized to the same maximum and 
are vertically displaced. The appearance and absence of beat- 
ing is discussed in the text. Compare with the corresponding 
theory (Fig.^J. 
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FIG. 9: (a) Measured optical THG intensity from sample B. 
Near T = 0, the THG spectrum clearly exhibits four peaks 
(the high-energy peak is amplified by a factor of 10). Normal- 
ized cuts at the spectral peak positions indicated by the white 
arrows in (a) are shown in (b) . Note the much slower, smooth 
decay of the envelope of the cut at 4.00 eV as compared to the 
cut at 5.16 eV, and the beating, which only occurs in the other 
two curves. 



tion peaks [see Fig.[5{b)] exhibit rather different spectral 
widths and our discussion of Sec. II VI The different spec- 
tral widths in Fig.EJa) correspond to strongly different 
decay times of the envelopes in Fig.^b). Again, only 
the envelopes of the first and fourth cuts show a smooth 
decay, whereas the envelopes of the second and third cuts 
exhibit a pronounced beating. 
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VI. CONCLUSIONS 

We have investigated the linear- and nonlinear-optical 
lineshapes of metal nanoparticles and metallic photonic 
crystal slabs. 

For particle ensembles, we have shown analytically and 
numerically that the comparison of time-resolved fem- 
tosecond second- or third-harmonic-generation experi- 
ments and extinction measurements does not allow one 
to distinguish between homogeneous and inhomogeneous 
contributions to the linewidth. Optical-rectification or 
four-wave-mixing experiments would provide such infor- 
mation. 

For metallic photonic crystal slabs, we have demon- 
strated that the model of two coupled Lorentz oscillators 
describes very well the key experimental features of lin- 
ear optics. In particular, Fano-like lineshapes appear in 
the absorption spectra. With regard to nonlinear op- 
tics, we have shown that - within the regime of "mod- 
erate coupling" - the nonlinear-optical third-harmonic- 
generation spectra provide information on the underlying 
source of the optical nonlinearity. Furthermore, the cal- 
culated nonlinear spectra reveal a beating in the spectral 
mixing products of the two peaks from linear optics, but 



not in the third harmonics of the latter peaks. 

Our corresponding experiments on third-harmonic 
generation of metallic photonic crystal slabs go beyond 
previous work regarding improved temporal resolution 
and the fact that we spectrally resolve the interferomet- 
ric third-harmonic signal. The spectra reveal a distinct 
behavior of the various spectral components versus time 
delay. Some spectral components exhibit a beating, oth- 
ers do not. Furthermore, the decay times of the envelopes 
strongly depend on the spectral component. The mea- 
sured spectra agree qualitatively very well with the pre- 
dictions of the simple theoretical model. The comparison 
allows us to identify the particle plasmon oscillation as 
the main source of nonlinearity. 

Acknowledgments 

We acknowledge support by the Center for Functional 
Nanostructures (CFN) of the Deutsche Forschungsge- 
meinschaft (DFG) within subproject A1.5. The research 
of M.W. is further supported by Project No. DFG- We 
1497/9-1. We thank H. Giessen, J. Kuhl, T. Zentgraf, 
and A. Christ for discussions. 



* Electronic address: Matthias.Klein@physik.uni-karlsruhe.de 

1 U. Kreibig and M. Vollmer, Optical Properties of Metal 
Clusters ( Springer- Verlag, Berlin, 1995). 

2 M. Neviere, in Electromagnetic Theory of Gratings, edited 
by R. Petit (Springer- Verlag, Berlin, 1980), Chap. 5, pp. 
123 - 157. 

3 S. Linden, J. Kuhl, and H. Giessen, Phys. Rev. Lett. 86, 
4688 (2001). 

4 A. Christ, S. G. Tikhodeev, N. A. Gippius, J. Kuhl, and H. 
Giessen, Phys. Rev. Lett. 91, 183901 (2003). 

5 T. Zentgraf, A. Christ, J. Kuhl, and H. Giessen, Phys. Rev. 
Lett. 93, 243901 (2004). 

6 B. Lamprecht, J. R. Krenn, A. Leitner, and F. R. 
Aussenegg, Appl. Phys. B: Lasers Opt. 69, 223 (1999). 

7 B. Lamprecht, G. Schider, R. T. Lechner, H. Ditlbacher, 
J. R. Krenn, A. Leitner, and F. R. Aussenegg, Phys. Rev. 
Lett. 84, 4721 (2000). 

8 B. Lamprecht, J. R. Krenn, A. Leitner, and F. R. 
Aussenegg, Phys. Rev. Lett. 83, 4421 (1999). 

9 G. Schider, J.R. Krenn, A. Hohenau, H. Ditlbacher, A. 
Leitner, F. R. Aussenegg, W. L. Schaich, I. Puscasu, B. 
Monacelli, and G. Boreman, Phys. Rev. B 68, 155427 
(2003). 

T. Vartanyan, M. Simon, and F. Trager, Appl. Phys. B: 
Lasers Opt. 68, 425 (1999). 

1 F. Stietz, J. Bosbach, T. Wenzel, T. Vartanyan, A. Gold- 
mann, and F. Trager, Phys. Rev. Lett. 84, 5644 (2000). 

2 M. J. Weida, S. Ogawa, H. Nagano, and H. Petek, J. Opt. 
Soc. Am. B 17, 1443 (2000). 

3 C. Sonnichsen, T. Franzl, T. Wilk, G. von Plessen, J. Feld- 
mann, O. Wilson, and P. Mulvaney, Phys. Rev. Lett. 88, 
077402 (2002). 

4 H. G. Craighead and G. A. Niklasson, Appl. Phys. Lett. 



44, 1134 (1984). 

15 M. Volmer and A. Weber, Z. Phys. Chem. (Leipzig) 119, 
277 (1926). 

16 B. Lamprecht, A. Leitner, and F. R. Aussenegg, Appl. 
Phys. B: Lasers Opt. 64, 269 (1997). 

17 The same problem associated with Eq. JIJ can arise for a 
conventional interferometric autocorrelation of laser pulses 
(Ref.iS), where one has P {2) (t) oc x (2) E 2 (t). There, 
however, the IAC remains unaffected, provided that the 
negative-frequency part of the laser spectrum has negligi- 
ble overlap with the positive- frequency part. 

18 Jean-Claude Diels and Wolfgang Rudolph, Ultrashort 
Laser Pulse Phenomena (Academic Press, London, 1996), 
Chap. 8. 

19 M. Wegener, D.S. Chemla, S. Schmitt-Rink, and W. 
Schiifer, Phys. Rev. A 42, 5675 (1990). 

20 A. Sharon, D. Rosenblatt, and A. A. Friesem, J. Opt. Soc. 
Am. A 14, 2985 (1997). 

21 U. Fano, Phys. Rev. 124, 1866 (1961). 

22 H. Giessen, S. Linden, A. Christ, J. Kuhl, D. Nau, T. 
Meier, P. Thomas, and S. W. Koch, in International Quan- 
tum Electronics Conference (IQEC), IFC5, San Francisco, 
2001 

23 Fitting the square root of the measured laser spec- 
trum (spectral intensity) with a sum of three Gaus- 
sians, we model the electric field spectrum with E(u) = 

ELi £ » e_( """° )2/ '^ where ^a/Si = 0.67, E 3 /Ex = 
0.65, ftoji = 1.36 eV, hu 2 = 1.67 eV, hu 3 = 1.82 eV, 
ha! = 0.11 eV, ha 2 = 0.17eV, and ha 3 = 0.048 eV (com- 
pare with gray area in Fig.|SJ. In the time domain, this 

corresponds to E(t) = ^ 3 =1 E„a„e~ (7 " t ^ cos(ui n t). 

24 A. Christ, T. Zentgraf, J. Kuhl, S. G. Tikhodeev, N. A. 



13 



Gippius, and H. Giessen, Phys. Rev. B 70, 125113 (2004). 
U. Morgner, F. X. Kartner, S. H. Cho, Y. Chen, H. A. Haus, 
J. G. Fujimoto, E. P. Ippen, V. Scheuer, G. Angelow, and 
T. Tschudi, Opt. Lett. 24, 411 (1999). 



M. U. Wehner, M. H. Ulm, and M. Wegener, Opt. Lett. 22, 
1455 (1997). 



